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Luttinger’s theorem is a fundamental result in the theory of interacting Fermi systems: it states that the volume 
inside the Fermi surface is left invariant by interactions, if the number of particles is held fixed. Although this 
is traditionally justified in terms of analytic properties of Green’s functions, it can be viewed as arising from 
a momentum balance argument that examines the response of the ground state to the insertion of a single flux 
quantum [M. Oshikawa, Phys. Rev. Lett. 84, 3370 (2000)]. This reveals that the Fermi volume is a topologically 
protected quantity, whose change requires a phase transition. However, this sheds no light on the stability or lack 
thereof of interacting semimetals, that either lack a Fermi surface, or have perfectly compensated electron and 
hole pockets and hence vanishing net Fermi volume. Here, I show that semimetallic phases in non-symmorphic 
crystals possess additional topological ‘Luttinger invariants’ that can be nonzero even though the Fermi volume 
vanishes. The existence of these invariants is linked to the inability of non-symmorphic crystals to host band 
insulating ground states except at special fillings. I exemplify the use of these new invariants by showing that 
they distinguish various classes of two- and three-dimensional semimetals. 


I. INTRODUCTION 

A well-trodden path to understanding many-electron sys¬ 
tems is to exploit proximity to a well-understood, typically 
free, model, and use this to access various properties of the 
interacting system through the lens of perturbation theory. 
There are few general statements made in this limit that re¬ 
main valid when the interactions are no longer parametrically 
small. A celebrated exception to this rule is Luttinger’s the¬ 
orem m it continues to impose constraints on the volume 
contained within the Fermi surface — the momentum-space 
surface that hosts low-energy excitations —even when the in¬ 
teractions are strong. Although Luttinger’s theorem is usually 
derived by examining the properties of Green’s functions to 
all orders in perturbation theory, about a decade ago Oshikawa 
gave it an elegant topological interpretation O : by determin¬ 
ing the change in symmetry of a many-electron ground state 
upon adiabatically inserting a single quantum of gauge flux, 
and computing the equivalent response of a Fermi liquid, the 
electron Ailing, a quantity that is determined by microscopic 
physics, can be related to the volume of the Fermi surface - a 
property of the low-energy effective theory. It is this linking 
of ultraviolet and infrared scales that makes Luttinger’s result 
a useful tool HID in analyzing correlated quantum matter. 

In recent years, the class of gapless Fermi systems in <7 > 1 
has grown to encompass semimetals such as graphene H) 
and its analogs in three dimensions EHID, where the Fermi 
energy is tuned to intersect a nodal surface where a pair 
of bands touch, either at discrete points or along continu¬ 
ous contours in the Brillouin zone. Understanding the prop¬ 
erties of such semimetals and their descendant phases is a 
major thrust of current research. The electronic dispersion 
near the nodal surface is typically linear, leading to a vanish¬ 
ing density of states at the Fermi energy, and consequently 
such systems are fairly robust: they are perturbatively sta¬ 
ble against interaction-induced band gaps, and any correlated 
gapped phases that emerge from these semimetals do so at 
fairly strong interactions. It is therefore desirable to place 
additional constraints on these systems, e.g. ones similar to 
those imposed by Luttinger’s theorem, that remain valid in 


the strong-coupling regime. However, semimetals have a van¬ 
ishing Fermi volume, and so Luttinger’s theorem is silent on 
their behavior. 

Two examples serve to illustrate why this is an unsatisfac¬ 
tory state of affairs. First, consider graphene’s honeycomb 
lattice at half-fllling. This corresponds to a single electron 
on each site, and therefore two electrons in each two-site unit 
cell. Famously, graphene is a semimetal at half-fllling: the 
Fermi surface consists solely of point nodes, and hence en¬ 
closes zero Fermi volume. An band theory analysis reveals 
that, absent interactions, it is impossible to open gap without 
breaking symmetry ifT^ . However, it is possible to construct 
a wavefunction that describes a gapped, symmetry-preserving 
phase of interacting electrons |[l2l[T3. This naturally raises 
the question of whether the existence of such a phase could be 
inferred on general grounds, rather than by laborious explicit 
construction — i.e., whether there is a Luttinger-like theo¬ 
rem for the nodal Fermi surface. A second example is pro¬ 
vided by compensated semimetals. These are systems whose 
Fermi surface consists of electron and hole pockets that en¬ 
close equal Fermi volumes, so that their net Fermi volume 
vanishes. Luttinger’s theorem, which is sensitive only to the 
net Fermi volume, is not readily applicable to this situation, 
that can arise in one of two ways. The first is to begin with 
a band insulator, and simultaneously pull the conduction (va¬ 
lence) bands below (above) the Fermi energy, while holding 
the particle number, and hence the Ailing, fixed. This ‘ac¬ 
cidental’ semimetal can be deformed back into an insulator 
by reversing this process. Another route begins with a nodal 
semimetal and pulls nodes above and below the Fermi energy. 
In this case, the system cannot be adiabatically deformed into 
a band insulator: its electron and hole pockets can at best be 
shrunk back to nodal points. How can we distinguish the two 
scenarios, sketched in Fig. simply by analyzing the Fermi 
surface and symmetries of a compensated semimetal? 

Here, I show that for certain semimetals, specifically those 
that occur in non-symmorphic crystals — loosely, crystals 
where screw axes or glide planes are essential to fully char¬ 
acterize the space group — one can associate additional topo¬ 
logical invariants with the Fermi surfaces (either nodal or 
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linger invariant. 

The remainder of this article is organized as follows. First, I 
introduce the flux insertion approach, and show how the non- 
symmorphic symmetry charge changes upon insertion of a 
single flux quantum. I illustrate the meaning of this change in 
terms of the spectral flow of a simple one-dimensional model 
of free fermions. I then turn to computing this spectral flow in 
a Fermi liquid, thereby linking the change in symmetry charge 
to a Fermi surface property. I then demonstrate how knowl¬ 
edge of the Luttinger invariant resolves the puzzles raised in 
the introduction, before closing with a discussion of applica¬ 
tions and extensions. 


II. FLUX INSERTION AND SPECTRAL ELOW 


FIG. 1. Two scenarios for compensated semimetals. Starting 
with equal-volume electron and hole pockets it may be possible to 
adiabatically deform the system into a band insulator (Case I), or 
there may be an obstruction to doing so, so that there remain one or 
more set of nodal points (nodal lines may also be possible) (Case II.) 
The Fermi surface in each case is shown inset. 


compensated). I demonstrate their existence via an analy¬ 
sis of crystalline point-group symmetry, focusing on the flow 
of discrete symmetry charge corresponding to a screw rota¬ 
tion or glide reflection into the system under a flux insertion. 
The result is an invariant that can be nonzero even when the 
Fermi sea volume vanishes. These invariants lead to rela¬ 
tions between low-energy parameters and high-energy prop¬ 
erties that are analogous to Luttinger’s theorem, and hence I 
dub them ‘Luttinger invariants’. Their existence is connected 
with a higher-dimensional analog of the Lieb-Schultz-Mattis 
theorem |[T4llT7ll that constrains the band structure of non- 
symmorphic crystals: except at certain speciflc integer flll- 
ings flxed by their space group such crystals can¬ 

not become insulating without either breaking symmetry or 
else exhibiting topological order, a phenomenon associated 
with an emergent gauge structure and the fractionalization of 
quantum. The flux insertion approach then reveals that the 
fractionalized quasiparticles that emerge in the topologically 
ordered phases transform non-trivially under the crystal sym¬ 
metry ED These arguments require that the fermions under 
consideration are spinless, or if spinful, possess an axis of spin 
conservation; they do not apply directly to spin-orbit coupled 
systems that lack this property. 

The Luttinger invariants so identifled can be used to dis¬ 
tinguish between different classes of gapless semimetals. For 
instance graphene possesses a vanishing Luttinger invariant, 
that can be computed directly from knowledge of symme¬ 
try. This immediately suggests that, with interactions, the 
semimetal can be gapped. Similarly, compensated semimet¬ 
als with (without) non-zero Luttinger invariants can (cannot) 
be deformed into a band insulator without breaking symmetry; 
we exemplify this dichotomy in d = 2 in a simple model of 
spinful electrons on the Shastry-Sutherland lattice, that hosts 
semimetallic phases both with and without a nonzero Lut- 


I begin by discussing the case of spinless fermions, before 
generalizing to spinful electrons (though without spin-orbit 
coupling). I consider a flnite crystal at Ailing o per unit cell, 
with primitive lattice vectors and periodic boundary con¬ 
ditions where Lai = 0, so that the system is deflned on a 
d-dimensional torus whose volume is unit cells. Note that 
here I specialize to the case of integer z^; as a result is an 
integer. The situation where o is not an integer was discussed 
previously O. Throughout, I work in units where h = e = 1, 
so that the quantum of flux, = 27r. 

The basic strategy is as follows: a single quantum of gauge 
flux is inserted adiabatically into the ground state of the sys¬ 
tem (alternatively, twisted boundary conditions are imposed 
around a non-contractible loop on the torus) and then a large 
gauge transformation is performed to to return the Hamilto¬ 
nian to its original form. Following Ref. El, I compute how 
symmetry quantum numbers of the ground state change in the 
resulting adiabatic cycle using two independent arguments. 
The first, ‘trivial symmetry counting’, is general and makes 
no assumptions as to the nature of the ground state — it de¬ 
pends only on the Ailing u. The second will make explicit 
reference to the ground state, by assuming it to be a Fermi 
liquid. Equating the results of these two ways of counting the 
change in symmetry charge, I arrive at a condition relating 
the Ailing to a Fermi surface parameter. It is this condition 
that depends crucially on the crystal structure. This argument 
closely parallels that of Oshikawa O, but differs from it in 
a key respect: rather than the change in momentum of the 
ground state, here more general symmetry operators lead to 
new topologically protected Luttinger invariants in the pres¬ 
ence of non-symmorphic symmetries. 

As a first step, I describe the flux insertion procedure. A 
necessary and sufficient condition so that the inserted flux is 
‘pure gauge’ for any non-contractible loop around the torus 
is to choose a reciprocal lattice vector K = rriibi, where 
the rrii are integers, and • bj = 27 T 6 ij, and then pick a 
gauge in which the Aharonov-Bohm flux ^ is rep¬ 
resented by the uniform vector potential A = ^ . I will 

assume further that the crystal is invariant under one or more 
non-symmorphic symmetry operations G = {g\T} that map 
points as G : r —> gr t, with gr = r. Note that a nec¬ 
essary condition for the operation to be non-symmorphic (TBi 
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is that the translation r is not in the lattice of discrete transla¬ 
tions, nor is it the projection of any lattice translation into the 
invariant subspace of g. Finally, I will assume that applying 
G consecutively Sq times is equivalent to the combination of 
a point-group operation and a lattice translation, i.e., is 
symmorphic; this defines the rank Sq of the operation G. 


A. Counting Symmetry Charge Microscopically 


For each non-symmorphic symmetry G = a specific 

fiux insertion that reveals the role of the symmetry can be con¬ 
structed as follows. Since G always involves a fractional lat¬ 
tice translation r, one can always choose a flux to insert such 
that the corresponding K is the smallest reciprocal lattice vec¬ 
tor parallel to r; note that it follows that K is left invariant by 
g, i.e., gK = K. Flux is then inserted adiabatically into 
the system by switching on a time-dependent vector potential 

~ ^ ~ 

end of the time T, there will be exactly one quantum of flux 
enclosed by a loop that encircles the torus parallel to K. The 
ground state |T^o) at t = 0 is assumed to preserve all the space 
group symmetries. Therefore it has a definite G-eigenvalue, 
given by G|T^o) = |T^o), where G is the unitary operator 

that implements the symmetry G on the Hilbert space. Since 
the Hamiltonian commutes with G at any t G [0,T] in this 
choice of gauge, it follows that at the end of the adiabatic flux 
insertion, the system is in some state |T^o) that has the same 
G-eigenvalue as the ground state, i.e. G|T^o) = e*^o|T^o). 
The final step is to return to the original gauge, by performing 
the large gauge transformation 


Uk = exp 




K • r)p{r) 


( 1 ) 


It is straightforward to demonstrate that 


UpGllK = Gexp 


T K 

I --- 



( 2 ) 


From this, it follows that the a full cycle of inserting a fiux and 
performing the gauge transformation ^ yields a state |^o) = 

I ^ 0)5 an eigenstate of G with eigenvalue where 


^iG ^ ^i{Go+uL^-^T-K) 


( 3 ) 


Eq. ^ gives the universal change in crystal symmetry charge 
upon a fiux insertion, independent of how this change is ac¬ 
commodated in the system, i.e., independent of the phase of 
the system. 


B. Spectral Flow: a one-dimensional example 

Before proceeding, a simple example may serve to illu¬ 
minate the meaning of the Luttinger invariant. Consider a 
one-dimensional lattice with two symmetries: translation (T^^) 
and a ‘non-symmorphic’ symmetry G that involves a discrete 
symmetry coupled with a half-translation — for instance, one 


could imagine a lattice with a two-site unit cell with reflection- 
conjugate orbitals on alternating sites (Fig. (e)). (Strictly 
speaking, there are no true non-symmorphic symmetries in 
d = 1, but this toy model serves to illustrate the general prin¬ 
ciple.) In units where the lattice spacing is a = 1, we have 
r = 2, and K = 27r. Now, consider non-interacting spinless 
fermions atu = 1 per unit cell, so that each site is half-filled. 

A glide mirror must square to a lattice translation: in the 
present example, G^ = Tx. The ID Bloch hamiltonian h{k) 
must therefore satisfy MG{k)h{k)MQ^{k) = h{k), where 
Mcik), which represents the action of the glide mirror on 
the Bloch states, is a function of k, a necessary condition for 
a glide. It is straightforward to see that Mcik) = 
where ax represents the action on the orbitals, and the 
ensures that [MG{k)]^ = = Tx in the Bloch basis. 

Thus, the Bloch states can be labeled by their eigenvalues 
m±{k) = Crucially, upon shifting k ^ k ^ 2 ti, 

the eigenvalues switch places, which requires that the two 
branches cross an odd number of times as the Brillouin zone is 
traversed (in higher dimensions, this would count the number 
of crossings along the glide plane.) With the further assump¬ 
tion that the model respects inversion symmetry, the cross¬ 
ing must occur at the zone boundary, k = tt ; on these gen¬ 
eral grounds, the band structure must take the form shown in 
Fig. (Note that the figure uses a slightly unconventional 
choice of the Brillouin zone, showing /c G [0, 27r] rather than 
the more familiar choice of [—7r,7r], to better illustrate the 
spectral flow.) At u = 1, the chemical potential lies at the 
crossing point of the two bands, so that the Fermi surface con¬ 
sists of a single point. 

Now, consider the spectral fiow induced by fiux inser¬ 
tion, which results in moving each single-particle state via 
k ^ k A At the nodal point (that we will denote ko 
in the general case), —k states from the valence band must 
evolve into Ak states in the conduction band, as fiux inser¬ 
tion preserves their G-eigenvalue. Similarly, —k states in the 
conduction band evolve into -\-k states in the valence band 
(Fig.j^a)). From this, it follows that the net result of the fiux 
insertion is to transfer a single charge from the conduction to 
the valence band (Fig.^b)) at the momentum k^^c = ^0 + x 
(In the thermodynamicHimit, /Cexc ^ ^o-)- Thus, after the fiux 
insertion we have = m_(/cexc)^+(^exc) = so 

that G — Go = TT (mod 27r), consistent with the microscopic 
counting in the preceding section for z/ = 1. Observe that a 
similar conclusion obtains as long as the bands cross an odd 
number of times. 

If it were possible to detach the bands without breaking 
symmetry (Fig. [^c)), then the corresponding spectral fiow 
would simply involve the fully filled valence band, and there 
would be no change in the symmetry quantum number of the 
ground state. The latter possibility is inconsistent with the 
microscopic counting: in other words, the non-symmorphic 
symmetry requires that the bands cross. Note however that 
at z/ = 2 , where both bands are filled, the adiabatic cycle 
is trivial and does not change the symmetry quantum num¬ 
ber (Fig. [^d)), but this is consistent with the counting from 
the microscopic theory. Thus, the microscopic counting for 
z/ = 1 constrains the band structure to require an odd number 
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FIG. 2. Spectral Flow at Integer Filling in c? = 1. (Brillouin zone 
shown in extended-zone scheme) (a) Free fermion dispersion for ID 
model with translation Tx and a single internal symmetry G involv¬ 
ing a half-translation. Bands are labeled by their G-eigenvalue m±i. 
Aiv — 1 all states are filled up to the nodal point, which is doubly- 
degenerate and therefore half-filled; upon flux insertion, filled quasi¬ 
particle states flow adiabatically (red arrows.), (b) Occupied states 
after flux insertion at z/ = 1: one charge has been transferred from 
valence to conduction band, leading to a shift in symmetry charge. 
Flux insertion leads to trivial spectral flow at z/ = 2 (c) or at z/ = 1 
when the symmetry G is broken and the bands are detached (d). 
A representative ID p-orbital lattice model with this symmetry is 
shown in (e) with the unit cell denoted by a red box. 


of crossings. 

This spectral flow picture may be generalized to d > 1 
by considering effective ID band structures for each distinct 
transverse momentum. We note that a similar discussion of 
the band structures required by non-symmorphic symmetries 
is given in Ref. 1^ . albeit without the discussion of spectral 
flow here. Ref. (241 discusses a closely related spectral flow 
on the surface of a topological crystalline insulator where one 
of the protecting symmetries is a glide mirror. 

An alternative way to phrase the above discussion is to sim¬ 
ply require that the the change in crystal symmetry quantum 
number computed in terms of the shift of quasiparticle states 
under /c —> ^ agrees with the microscopic counting. If the 
microscopic counting involves a nontrivial symmetry change, 
then it follows (by arguments similar to the example above) 
that this must be reflected in the low-energy description. If 
the system is in a Fermi liquid phase, this gives an invariant 
associated with the Fermi surface response to flux insertion. 
I now adopt this approach to compute a ‘Luttinger invariant’ 
associated with a Fermi surface of zero volume for (i > 1. 


C. Counting Symmetry Charge in the Fermi Liquid 

I now compute the flow of symmetry charge upon flux in¬ 
sertion assuming the system is in a Fermi liquid phase, and 
so can be described in terms of long-lived quasiparticles, with 
low-energy effective Hamiltonian 

ZZ ~ y] e(fe)nfe + /(fe, , (4) 

k k,k' 

Assume that the ground state before flux insertion is the fllled 
Fermi sea of quasiparticle states. Under the flux insertion, 
any change in the ground state can be accounted for in terms 
of quasiparticle excitations generated near the Fermi surface. 
Recall that quasiparticles near the Fermi surface are long- 
lived; using this along with the adiabatic continuity with the 
Fermi gas, one may determine the change in the quasiparticle 
population Shk. Observe that flux threading for noninteract¬ 
ing fermions simply shifts the Fermi sea by a uniform amount 
^, as each k ^ k Now, a uniform shift of the Fermi 

surface of this form amounts to producing quasiparticles on 
one side of the Fermi surface, and quasiholes on the opposite 
side. Therefore, as the excitations are close to the Fermi sur¬ 
face, Fermi liquid theory may be used to compute the change 
in the ground-state properties. 

For this purpose, it is convenient to split the behavior of 
the glide operation into its two distinct components: the 
point group operation and the fractional translation, by writ¬ 
ing G = gT^-. Consider the states iT^o)^ l^o) before and after 
flux insertion. Then, (using similar notation to the previous 
section) 


= (folflT^l^o). (5) 

Now, observe that the shift of the Fermi surface was by a 
^-invariant momentum (recall gK = K), so it follows that 
the only difference between G and Go emerges from the frac¬ 
tional translation by r. Since the result of the adiabatic pro¬ 
cess is the shift of the whole Fermi sea by K/L, the change in 
the momentum P of the system during the adiabatic process 
is given by (2]-[4l 


ap = n^Ay 

where is the integer ‘occupation’ number of fllled quasi¬ 
particle states in the ground state. (I note that a more pre¬ 
cise computation actually obtains AP by first summing over 
deformations of the quasiparticle distribution near the Fermi 
energy, converting the sum to an integral over the Fermi sur¬ 
face, and then using Stoke’s theorem to relate this to the Fermi 
surface volume, and hence the number of fllled quasiparticle 
states relative to the bottom of the conduction band lHHH. As 
this yields the same result for L ^ oc, I simply use the sim¬ 
ple formula above.) From and the preceding discussion 
it is possible to compute the change in symmetry charge in a 
Fermi liquid: 


^ i(Go+ivW^) 


(7) 
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where we have used the fact that the fractional translation has 
eigenvalues An explicit calculation of the symmetry 

charge in the free fermion case is provided in an appendix. 

Eq. (|^ expresses the change in the symmetry charge assum¬ 
ing that the low-lying excitations form a Fermi liquid. Com¬ 
paring this to the microscopic counting in ([^, yields a consis¬ 
tency condition, namely that ■ K is 

an integer multiple of 2tt. 

To proceed, one must determine the value of r • K. This 
depends crucially on the fact that G is non-symmorphic. From 
the definition of the rank, is symmorphic, so that the 
associated translation Sqt must be a lattice translation. Since 
K is the smallest reciprocal lattice vector parallel to K, it 
follows at that t K = 27tpISg foi* some integer p < So, 
relatively prime to So- From this, 

p{N^P - vL"^) = Lx (integer) x Sg- (8) 

Let us take L relatively prime to Sg, and define ^ = 
so that 

p(Xf ^ = L X (integer) x Sq- (9) 

Since the number of filled quasiparticle states in the system is 
extensive, it follows that as L ^ oc, is independent of L. 
Since u is also independent of L, (|^ cannot be consistent in 
the thermodynamic limit unless the RHS also scales with 
L^. Thus, canceling a factor of from both sides of ([^ 
leads to p{x^p^ ~ — (integer) x Sq- Finally, as p and Sg 

are relatively prime, the only way to satisfy this relation is if 

Xf = — nSc (10) 

for n any integer (I have dropped the superscript on with 
the understanding that I refer to its L ^ oc limit.) 

III. LUTTINGER INVARIANT 

In a crystal whose space group contains many non- 
symmorphic symmetries, one may obtain a similar constraint 
on xf for each such symmetry. From this and the fact that n 
can be any integer, it follows that if 5* is the the least common 
multiple of all the Sg, then 

Xi? = (mod 5*) (11) 

is a topological invariant of the system, that I propose to 
term a ‘Futtinger invariant’. For any crystal containing non- 
symmorphic operations, 5* > 1; of the 157 non-symmorphic 
space groups, 155 are of this type. However, there exist two 
‘exceptional’ non-symmorphic space groups where every in¬ 
dividual operation can be rendered symmorphic by a change 
in real-space origin, so that the argument above does not ap¬ 
ply immediately. Still, for these groups it is possible to show 
suitable combinations of consecutive fiux insertions yield sim¬ 
ilar constraints, so that they too have 5* > 1 (both have 
iS* = 2). Therefore, 5* > 1 if and only if the crystal is 
non- symmorphic. 


Eq. ([n) (or equivalently, ( p^ ) is the central result of this 
paper: it relates the filling, defined in the microscopic theory, 
to the counting of quasiparticle excitations in the emergent 
low-energy description Q. The second term on the right-hand 
side of ( p^ refers to the filled bands; in essence, this equation 
refiects the fact that all the charge in the system is either bound 
up into filled bands or contributes to the low-energy gapless 
quasiparticle excitations. The presence of 5* in this expres¬ 
sion is linked to a topological requirement that energy bands 
in crystals can only appear in multiplets (containing a multi¬ 
ple of 5* bands) that ‘stick’ together f\M . Note that there is 
a subtle distinction between the non-symmorphic rank 5 of a 
space group as defined in Ref. and the definition of 5* 
used here; I comment on this in the appendix, but this distinc¬ 
tion is unimportant for the examples studied here. 

An expression similar to Eq. ^TT\ was derived in Ref. Q: 
if Vf = {2fYxf, denotes the volume of the Fermi sea, then 
— V must be an integer in a translationally-invariant sys¬ 
tem. For fractional v, this is the familiar result that the Fermi 
sea volume is protected: it is proportional to the fractional part 
of the filling. When the filling is an integer (say f = 1, for 
specificity) the volume of the Fermi sea vanishes. This is in 
accord with the intuition (to be qualified shortly) that a sys¬ 
tem with all states in the Brillouin zone filled is inert and can 
be adiabatically deformed into one with no gapless excitations 
and hence no Fermi surface. In this manner, Futtinger’s the¬ 
orem can be seen to have a topological origin: the volume of 
the Fermi sea is an invariant that is ‘protected’ against inter¬ 
actions as long as the low-energy description takes the form 
of Eq. 0 . Extending this picture to include more general 
symmetries beyond translation, it is clear that 5* = 1 , corre¬ 
sponding to a symmorphic crystal, is uninteresting: there are 
no additional Futtinger invariants besides the Fermi sea vol¬ 
ume. When this vanishes, there is indeed no obstruction to 
adiabatically deforming the system into a gapped phase. 

Non-symmorphic crystals where 5* > 1 possess additional 
invariants that modify this picture. From ^TT) one can see 
that as long as u is indivisible by 5*, it must be true that 
Xf = z^(modtS*) 7 ^ 0. As we have noted, for any inte¬ 
ger F, the volume of the Fermi sea vanishes in the Brillouin 
zone, so xf must be a new invariant distinct from the Fermi 
sea volume. Observe that xf counts the symmetry change 
induced by low-energy quasiparticles described by Q; thus, 
the fact that Xf is non-zero means that these excitations must 
be gapless. (For, if they were gapped, it would be possible 
to smoothly change parameters until they were very far away 
from the Fermi surface and therefore could not contribute to 
Xf-) One way to reconcile the vanishing of the Fermi sea vol¬ 
ume with the gaplessness of Q is if the Fermi surface con¬ 
sists entirely of point or line nodes, i.e., describes a nodal 
semimetal. An alternative is a compensated system, with 
equal-volume electron and hole Fermi surfaces, although this 
may be ruled out by additional symmetries ll^ . Rather than 
the momentum balance that leads to the protection of a Fermi 
sea, here the topological protection is linked to the quantum 
numbers of discrete spatial symmetries such as rotation or re- 
fiection. Thus, non-symmorphic symmetries allow us to iden¬ 
tify a class of semimetals that are protected by a topological 



6 


‘Luttinger’ invariant analogously to how a filled Fermi sea is 
protected by Luttinger’s theorem. 

Note that, as in the case with a nonzero Fermi sea volume, 
one way to evade the protection is that the low-energy the¬ 
ory no longer takes the form of a Fermi liquid so that Q is 
no longer a valid description of the system. However, even a 
non-Fermi liquid without sharply-defined quasiparticles, that 
nevertheless has an appropriately defined Fermi surface of 
low-energy excitations, is still amenable to this argument and 
therefore can have a nonzero Luttinger invariant; similar re¬ 
sults then apply. There is a familiar example in d = 1: Lut¬ 
tinger liquids satisfy Luttinger’s theorem, but are not Fermi 
liquids (271 [28l. Alternatively, the system could open a gap, 
but cannot enter a trivial insulating phase (meaning one whose 
wave function can be adiabatically continued to that of a band 
insulator) without breaking symmetry fTsl , refiecting the gen¬ 
eralized Lieb-Schultz-Mattis theorem. Note that all these ar¬ 
guments require the presence of f/(l) charge conservation and 
the space group symmetries; however, the breaking of these 
symmetries is usually detectable in experiments. 

This picture shows the crucial role played by the non- 
symmorphic symmetry: on the one hand, it has a nontrivial 
relationship with translation, meaning that adiabatic shifts of 
momenta change its value even for integer u; on the other 
hand, it is defined along a line or plane in the Brillouin zone, 
and therefore we can track the spectral fiow of quasiparti¬ 
cles in a concrete way by assigning definite symmetry quan¬ 
tum numbers to bands along these high-symmetry directions. 
Note that the analysis of band structures is limited to non¬ 
interacting systems, but the arguments that led to the compu¬ 
tation of the Luttinger invariant are non-perturbative in nature, 
and so apply more generally. The non-perturbative Luttinger 
invariant ([n) provides a simple way to compute, purely from 
knowledge of stoichiometry and crystal structure, whether a 
given semimetallic dispersion must exhibit a nontrivial spec¬ 
tral fiow even in the presence of interactions, and hence iden¬ 
tify protected semimetals. 


A. Including Spin 

The extension of these arguments to electrons with spin is 
straightforward in the absence of spin-orbit coupling. Since 
the two spin species are independently conserved, one may 
define separate fillings so that u = Accordingly, it is 

possible to introduce a fictitious gauge flux that couples to the 
up and down spins separately, and follow the reasoning above 
to obtain a pair of invariants, Xf — (mod 5*) with a =t 
, In the spin symmetric case where we may simply 

write Xf = (mod 5*): even though there are 

twice as many electrons in each unit cell because of the spin 
degeneracy, in this case the topological invariant is computed 
by simply halving the filling and therefore coincides with that 
computed for spinless fermions. Note that, in applying this 
argument, we assume that the spin transforms trivially under 
the mirror component of the glide symmetry — a choice that 
is only consistent absent spin-orbit coupling, that breaks the 
spin rotation symmetry. 


IV. EXAMPLES 

We now revisit the examples from the introduction, armed 
with the new Luttinger invariant. First, consider graphene, 
where the spin-orbit coupling is negligible and can be set to 
zero. At half-filling on the honeycomb lattice, v = = 

2, owing to the two sites in each unit cell; therefore, the Fermi 
volume vanishes. As the honeycomb lattice has a symmorphic 
space group {P^lmmc), it has no other nontrivial Luttinger 
invariants; and hence no generalized ‘Luttinger theorem’ can 
be associated with Dirac nodes in graphene. This is consistent 
with the existence of a gapped, symmetry-preserving phase of 
spinful fermions at f = 2 on the honeycomb lattice with no 
fractionalization (121 [TSl ; such a phase cannot descend from 
a gapless system with a nonzero Luttinger invariant. Thus, 
knowledge that the Luttinger invariant vanishes allows one to 
deduce the existence of a gapped symmetry-preserving phase 
without explicit construction of a microscopic model or trial 
wavefunction. 

Turning to our second puzzle of compensated semimetals, it 
is useful to focus on an illustrative example in d = 2 provided 
by electrons on the Shastry-Sutherland lattice (SSL) (29l . 
again without spin-orbit coupling. The SSL has a non- 
symmorphic space group ip^g) with 5* = 5 = 2. Us¬ 
ing a simple s-orbital tight-binding model we find that the 
non-interacting band structure may exhibit nodal semimetal¬ 
lic behavior at even integer values of f = f>^ + f^. Con¬ 
sider the half-filled case, where f = 4, so that f^ = fj^ = 2 
(Fig.j^b)). From the preceding arguments, at this filling there 
is no Luttinger invariant associated with the energy bands: 
Xf = (mod 5*) = 0. Indeed, it is straightforward to write 
down a simple insulating wavefunction at this filling ( 23 , by 
binding electrons into singlets placed on the solid bonds in 
Fig. [^a); alternatively, suitable choices of the hopping pa¬ 
rameters leads to a gap opening (Fig. [^c)) without a sym¬ 
metry change. On the other hand, at quarter filling (f = 2, 
~ considerations presented here reveal that 

the semimetallic phases are protected by a non-zero Luttinger 
invariant (xf = 1), and therefore cannot be gapped with¬ 
out triggering fractionalization or breaking symmetry. While 
semimetallic phases that appear in the SSL have been studied 
recently (30l, the essential distinction between those at f = 2 
and F = 4 seems to have been overlooked. 

Similar statements can be made for three-dimensional 
non-symmorphic crystals, e.g^ the diamond or hexagonal 
close-packed structures. Specifically, any filling-enforced 
semimetal ED in a non-symmorphic crystal is associated with 
a non-zero Luttinger invariant. To see this, we observe that a 
filling-enforced semimetal occurs in a non-symmorphic crys¬ 
tal at any filling (of spinful electrons) that is not an even inte¬ 
ger multiple of the non-symmorphic rank; the non-zero value 
of the Luttinger invariant then follows immediately from the 
arguments of the preceding section. Technically, our iden¬ 
tification of the Luttinger invariant relies on the presence of 
spin-rotation symmetry and hence the absence of spin-orbit 
coupling, that may be present for a generic filling-enforced 
semimetal. However, motivated by the absence of symmetry¬ 
preserving, non-fractionalized insulating ground states at such 


7 


(a) ^ ^ ^ (^) 2t' > t (c) 2t' < t 




Xf = 1 Xf = 0 Xf = 1 Xf = 1 Xf 1 

FIG. 3. Luttinger Invariants in the Shastry-Sutherland Lattice. 

(a) Shastry-Sutherland lattice with pAg space group and Brillouin 
zone with high-symmetry directions labeled. (b,c): Tight-binding 
band structure and Fermi surfaces at even integer fillings for repre¬ 
sentative values of t'/t\ electron (hole) pockets are colored in red 
(blue). Owing to the combination of inversion and time-reversal 
symmetry, each energy level is two-fold degenerate. For 2t' > t, (b) 
the system is a (perfectly compensated) semimetal at one-quarter-, 
half- and three-quarter- filling {u = 2,4, 6 for spinful electrons.). 
For V — 2,6 the Fermi surface consists of electron and hole pockets 
enclosing zero net area, whereas z/ = 4 it consists of a quadratic band 
crossing at the zone center and a non-dispersing line of filled elec¬ 
tronic states connecting the zone center to the corners. The semimet¬ 
als at z/ = 2,6 semimetals are protected by a nonzero Luttinger in¬ 
variant (xf = 1) and remain gapless as we tune parameters with¬ 
out changing symmetry, but the unprotected nodes can be gapped, 
as show] in (c), where 2t' < t. In the strongly interacting limit 
at half-filling the effective description is a Heisenberg model whose 
(unique) ground state is a crystal of singlets on strong bonds 
(solid lines in (a)); this phase can be adiabatically connected to the 
z/ = 4 band insulator in (c). the non-zero Luttinger invariant guar¬ 
antees that any gapped symmetric ground state at z/ = 2,6 has topo¬ 
logical order. Note that the fact that the electron and hole pockets 
themselves touch in the Brillouin zone is an artefact of the neglect of 
spin-orbit and higher-neighbor couplings. 



fillings even with spin-orbit interactions ||20l, I conjecture that 
the Luttinger invariant remains non-zero also in this case. 


V. CONCLUDING REMARKS 


In closing, I comment briefly on the applicability of these 
ideas to three dimensional nodal semimetallic phases. Re¬ 
cently, Dirac semimetallic phases have been identifled in non- 
symmorphic crystals in two 1^ and three (T] [31] [32]| di¬ 
mensions. The corresponding band structures would have 
non-zero Luttinger invariants in the limit of vanishing spin- 
orbit coupling; however, in this limit the bands can touch 
along nodal surfaces rather than at isolated points, and it is 
not clear whether the invariants as deflned here survive the 
inclusion of spin-orbit coupling. Note that a quite different 
approach 1 ^ than flux-threading seems necessary to extend 
similar arguments for gapped phases to systems without spin- 
rotation symmetry; whether such techniques can be suitably 
adapted to treat gapless systems remains an outstanding prob¬ 
lem. Therefore, whether a similar invariant can be identifled 
in the presence of strong spin-orbit coupling is at present an 


open question, that seems worthy of further study. 

It is possible to extend similar ideas to Kondo lattice mod¬ 
els believed to capture the essential physics of heavy-fermion 
materials. Here, a generalized Luttinger invariant can be com¬ 
puted if the Ailing is taken to count both the conduction elec¬ 
trons and the local spin moments within the unit cell. The 
derivation of this invariant and its consequences for heavy 
fermion systems are beyond the scope of the present article, 
but will be discussed elsewhere 1 ^1^ . 

Finally, I emphasize that while Dirac semimetals with 
nodes on their zone boundaries do not generically have 
topological surface states or a nontrivial bulk electromagnetic 
response — and are thus not ‘topological’ in one sense that 
is currently in vogue — the presence of a nonzero Luttinger 
invariant means that they are parent semimetals for three- 
dimensional topologically ordered (i.e., fractionalized) phases 
that emerge when they are gapped while preserving symme¬ 
try. This suggests that understanding such semimetals and 
their instabilities is one possible route to new phases of matter. 
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Appendix A: Counting Symmetry Charge for Free Fermions 

Although the argument in the main body of the paper that 
restricted consideration to the low-energy excitations near the 
Fermi surface (where there are well-deflned, scattering free 
quasiparticles), the conclusion was that the change in symme¬ 
try charge upon flux insertion was simply related to the num¬ 
ber of filled momentum-space states. One might intuitively 
(though non-rigorously) argue that this quantity could be di¬ 
rectly computed for free fermions, and that its value is an adi¬ 
abatic invariant as interactions are switched on. 

Let us determine the change in the symmetry charge of the 
system in terms of occupied (single-particle) fermionic states 
in the free Fermi gas. First, note that the quasiparticle creation 
operator transforms under the symmetry G = {^flr} as 

(Al) 

where Before inserting flux, the system contains 

a fllled Fermi sea: |T^o) = life occ show 

that 


G\^o)= n (A2) 

k,occ. 

where I have used the fact that the set of all fllled states must 
be invariant under G, for otherwise the fllled Fermi sea would 
break symmetry. Upon inserting the flux, the momentum of 
each quasiparticle state has shifted sls k —> /c + so that 
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l^o) = rife.occ. 4+^ therefore, 


G|^o) = e*^'«,oc=.('^-'“+^)|#o) = e*'^|^o)- (A3) 


If one defines Go = occ ^ ‘ from (A3), 

giG ^gi(Go+iV^^>^)^ 

which reduces to the result 0 obtained in the main text. 


(A4) 


Appendix B: Relation between S* and Non-Symmorphic Rank 

As noted in the main text, there is a subtle distinction be¬ 
tween the non-symmorphic rank 5 of a space group as defined 
in Ref. ifT^ and the definition of 5* used here. While these 


coincide in most cases, the strict definition of the rank is the 
minimum filling v at which a symmetric insulating state is 
possible without triggering fractionalization, and is hence not 
obviously related to the definition of 5*. In spite of this sub¬ 
tlety, the fact that bands always appear in multiplets of S in 
non-symmorphic crystals, suggests that by the heuristic argu¬ 
ment on accommodating charge in filled bands versus gapless 
quasiparticles, S should replace 5* in the above expressions. 
However, the precise definition of the rank requires one to 
consider the electronic polarization, a quantity that is well- 
defined for an insulator but not in the present case, where 
there are gapless excitations that carry charge. Therefore, I 
am unable at present to give a rigorous justification to replace 
5* by S in ((n), and can only conclude that S must at least 
be divisible by 5* and that 5* > 1 if and only if the crystal 
is non-symmorphic. In practice, however, a plurality of non- 
symmorphic crystals have 5 = 2, and therefore 5* = 5 in 
these cases. 
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